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Abstract 

We study the error rate of CNOT operations in the Kane solid state 
quantum computer architecture pQ. A spin Hamiltonian is used to de- 
scribe the system. Dephasing is included as exponential decay of the 
off diagonal elements of the system's density matrix. Using available spin 
echo decay data, the CNOT error rate is estimated at approximately 10~ 3 . 

1 Introduction 

Existing classical computers manipulate bits that can be exclusively or 1. 
Quantum computers manipulate two level quantum systems called qubits that 
can be arbitrary superpositions of |0) and While the idea of a quantum 
computer was first suggested by Benioff and Feynman in the early 80s 012], 
the first quantum algorithm that could solve an interesting real world problem 
faster than it's classical equivalent was published by Shor in 1994 4 . Shor's 
algorithm factorizes integers with the number of steps growing polynomially in 
the number of digits whereas the best known classical algorithm grows expo- 
nentially. A significant milestone in the quest to build a quantum computer was 
the construction of a 7 qubit liquid NMR implementation of Shor's algorithm 
designed to factorize 15 Unfortunately, the liquid NMR approach is not 
expected to work beyond a few tens of qubits 

Many different technologies are being researched in the hope of producing a 
scalable quantum computer. A sample of the diverse proposals can be found in 
PEUHlinilinilll]- In Kane's solid state proposal [DIH], the nuclear spins of 
single 31 P dopant atoms in 28 Si are used as qubits. This approach aims to take 
maximum advantage of the industry expertise acquired during the last 50 years 
of conventional semiconductor electronics. 

In this paper we study the error rate of CNOT operations in the Kane 
quantum computer. Initial simulations were carried out without dephasing to 
enable the pulse profiles of the controlling electrodes to be optimized ^3]. The 
lowest error rate achieved in the absence of dephasing was 5 x 10~ 5 . When a 
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physically reasonable level of dephasing was included in the simulation, the error 
rate increased to approximately 10~ 3 . While theoretical estimates of the error 
rate required for fault tolerant computation are of order 10 -6 |14|. numerical 
simulations by Zalka suggest an error rate of 10~ 3 or higher may be tolerable 
15 . Further work is required to determine the maximum allowed error rate in 
the Kane architecture. 

This paper is organized as follows. In section [2] the physical architecture of 
the Kane quantum computer is described. In section|21the process of performing 
a CNOT operation in the Kane architecture is presented with emphasis on 
achieving the lowest possible error rate in the absence of dephasing. Further 
details of this process can be found in In section 0] our model of dephasing 
is described and its effect on the error rate of the CNOT gate. In section [5] 
we conclude with a discussion of the implications of our estimate of the likely 
minimum error rate of the Kane CNOT gate. 

2 The Kane Quantum Computer 

The 31 P in 28 Si system is thought to be well suited to use as a qubit due to it's 
long relaxation (Ti) and dephasing (T2) times. Both times only have meaning 
when the system is in a steady magnetic field. Assuming the field is parallel with 
the z-axis, the relaxation time refers to the time taken for 1/e of the spins in 
the sample to spontaneously flip whereas the dephasing time refers to the time 
taken for the x and y components of a single spin to decay by a factor of 1/e. In 
natural silicon containing 4.7% 29 Si, relaxation times in excess 1 hour have 
been observed for the donor electron at T=1.25K and B^0.3T |17|. The nuclear 
relaxation time has been estimated at over 80 hours in similar conditions |18|. 
The donor electron dephasing time T2 in enriched 28 Si containing (0.12±0.08)% 
29 Si mi] has been measured at T=1.4K to be ~0.5ms [20] - At the time of writing, 
no experimental data relating to the nuclear dephasing time has been obtained 
to the authors' knowledge. 

The phosphorous donor electrons are used primarily to mediate interactions 
between neighboring nuclear qubits. As such, they are polarized to remove their 
spin degree of freedom from the system. This can be achieved by maintaining a 
steady B Z =2T at around T=4K ^B]. To take advantage of the long Ti and T 2 
times discussed in the previous paragraph, the operating temperature will more 
likely need to be ~1K. Techniques for relaxing the high field and low temperature 
requirements such as spin refrigeration are under investigation |12) . 

In the Kane architecture, qubits are arranged in a single line. Control is 
achieved via electrodes above and between each qubit and a global transverse os- 
cillating field of magnitude ~10~ 3 T (Fig.QJ. To selectively manipulate a single 
qubit, the A-electrode above it is biased. A positive/negative bias draws/drives 
the donor electron away from the nucleus reducing the magnitude of the inter- 
action between the electron and nuclear spins. This in turn reduces the energy 
difference between nuclear spin up (|0)) and down (|1)) allowing this transition 
to be brought into resonance with a globally applied oscillating magnetic field. 
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B ac ~10- 3 T 



Figure 1: Schematic of the Kane architecture. The rightmost 2 qubits show the 
notation to be used when discussing the CNOT operation. 



Depending on the timing of the A-electrode bias, the qubit can be rotated into 
an arbitrary superposition a|0) + /3|1). Clearly this scheme also allows arbi- 
trary combinations of individual qubits to be simultaneously and independently 
manipulated. 

Interactions between neighboring qubits are governed by the J-electrodes. 
A positive bias encourages greater overlap of the donor electron wave functions 
leading to indirect coupling of their associated nuclei. In analogy to the single 
qubit case, this allows two qubit transitions to be performed selectively between 
arbitrary neighbors. A discussion of the electrode pulses required to implement 
a CNOT gate is given in the next section. 



3 The CNOT Gate on a Kane QC 

Performing a CNOT operation on a Kane QC is an involved process described in 
detail in [T5|. Given the high field (2T) and low temperature (~1K) operating 
conditions, we can model the behavior of the system with a spin Hamiltonian. 
Only two qubits are required to perform a CNOT operations so for the remainder 
of the paper we will restrict our attention to a computer with just two qubits. 
The basic notation is shown in (Fig. Furthermore, let cr^ = a z ®I®I®I, 
cr z el = I ® a z ® I ® I, a z n2 = I ® I ®a z ® I and a z e2 = I®I®I®a z where / is 
the 2x2 identity matrix, a z is the usual Pauli matrix and ® denotes the matrix 
outer product. With these definitions the meaning of terms such as er^ 2 and o e \ 
should be self evident. 

Let g n be the g-factor for the phosphorus nucleus, [i n the nuclear magneton 
and /ib the Bohr magneton. The Hamiltonian can be broken into three parts 

H = H Z + H int (t)+H ajc (t). (1) 
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The Zeeman interactions terms are contained in Hz 

Hz = -g n ^ n B z {a z nl + a z n2 ) + fi B B z (a z el + a z e2 ). (2) 

The contact hyperfine and exchange interaction terms, both of which can be 
modified via the electrode potentials are 

-ffint(i) = Ai(t)a n i ■ o ei + A 2 {t)3 n2 ■ a e2 + J(t)a e i ■ a e2 , (3) 

where Ai(t) = 87r^sg„^„|$ i (0)| 2 /3, |3>i(0)| is the magnitude of the wavefunc- 
tion of donor electron i at phosphorous nucleus i, and J(t) depends on the 
overlap of the two donor electron wave functions. The dependance of these 
quantities on their associated electrode voltages is a subject of ongoing research 
ESI- In this paper the hyperfine and exchange interaction magnitudes 
Ai and J will frequently be discussed as though directly manipulablc. 

The last part of the Hamiltonian contains the coupling to the global oscil- 
lating field B ac - 

ffac(i) = B ac (t) cos(w«) + a x n2 ) + [i B Wi\ + o*n)\ 

+ B„(t) sm(ut)[-g nf i n (a y nl + a v n2 ) + fi B (a v el + a v e2 )]. (4) 

Using the above definitions, only the quantities At, J and -B ac need to be ma- 
nipulated to perform a CNOT operation. 

For clarity assume the computer is initially in one of the states |00), |01), 
1 10) or 1 11) and that we wish to perform a CNOT operation with qubit 1 as the 
control. Step one is to break the degeneracy of the two qubits' energy levels to 
allow the control and target qubits to be distinguished. To make qubit 1 the 
control the value of Ai is increased (qubit 1 will be assumed to be the control 
qubit for the remainder of the paper). 

Step two is to gradually apply a positive potential to the J electrode in order 
to force greater overlap of the donor electron wave functions and hence greater 
(indirect) coupling of the underlying nuclear qubits. The rate of this change is 
limited so as to be adiabatic — qubits initially in energy eigenstates remain in 
energy eigenstates throughout this step. 

Let |symm) and |anti) denote the standard symmetric and antisymmetric 
superpositions of 1 10) and |01). Step three is to adiabatically reduce the A\ 
coupling back to its initial value once more. During this step, anti-level-crossing 
behavior changes the input states as 1 10) — > |symm) and |01) — > |anti). 

Step four is the application of an oscillating field -B ac resonant with the 
|symm) <-» 1 1 1 ) transition. This oscillating field is maintained until these two 
states have been interchanged. Steps five to seven are the time reverse of steps 
one to three. The process is shown schematically in (Fig. El- Note that steps 1 
and 7 (the increasing and decreasing of A^ ) have been omitted as the only limit 
to their speed is that they be done in a time much greater than ?i/0.01eV~ O.lps 
where O.OleV is the orbital excitation energy of the donor electron. 

In general the fidelity of the adiabatic steps in the procedure can be increased 
arbitrarily by making them indefinitely long. In reality, of course, this is not 
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Figure 2: Gate profiles and state energies during a CNOT operation in units of 
g n \x n B z = 7.1 X 10- 5 meV. 
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practical and the procedure must be implemented on a time-scale that is short 
compared to the systems dephasing time. For a given H(t), the degree to which 
the evolution deviates from perfect adiabaticity can quantified by the measure 



0(t) = Max a#b 



h\{Mt)\&mt))\Mt))\ 



({Mt)\H(t)\Mt)) - mt)\H(t)\MtW 



(5) 



It is desired that 0(i) <C 1. Here the states \ip a (t)) are the set of eigenstates of 
the Hamiltonian at time t. It is possible to reduce 0(t) without increasing the 
duration of the step by optimizing the profiles of the time dependent parameters 
in the Hamiltonian. In the case of the Kane architecture this means optimizing 
the shape of the evolution of Ai(t) and J(t). 




Figure 3: Possible forms of the J(t) profile for step 2 of the adiabatic CNOT 
gate. J(t) is in units of g n ^nB z = 7.1 x 10~ 5 meV. 

Various profiles for the adiabatic steps in the CNOT procedure have been 
investigated in ^3]. In (Fig. [3}, we have plotted three possible J(t) profiles for 
step 2 of the CNOT gate. The function O(i) for each profile is shown in (Fig.^J. 
Profile 1 is a simple linear pulse, profile 2 can be seen to be the best of the three 
and is described by J(t) — 810a(l — sech(5i/r)) where r = 9/is is the duration 
of the pulse and a — 1.0366 is a factor introduced to ensure that J(r) = 810. 
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Figure 4: The adiabatic measure Q(t) for each J(i) profile. 



The third profile 
J(t) = I ^ 



t(l+7r/2) 



T ' 

1 + sin 



7T t-T/(l+7T/2) 

2 T/(1+2/tt) 



0<t <T/(1 + tt/2) 
T/(1 + tt/2) < t < T 



(6) 



although not quite as efficient as profile 2, is a composite linear-sinusoidal profile 
that was used in the calculations presented in this paper due to numerical dif- 
ficulties in solving the Schrodinger equation for profile 2. The advantage of the 
second two profiles over the linear one is that they flatten out as J approaches 
810. At J = 816.65, the system undergoes a level crossing. To maintain adi- 
abatic evolution, J(t) needs to change more slowly near this value. Note that 
the reason it is desirable to make J(t) so large is to ensure that there is a large 
energy difference between |symm) and |anti) during step 4 (the application of 
B ac ). This difference is given by 

5E = 2A 2 (- ' 



f 



(7) 

M-B-Bz + 9n^ n B z ^bB z + g n [i n B z - 2 J 

Without a large energy difference, the oscillating field B, AC which is set to res- 
onate with the transition |symm) <-> |11) will also be very close to resonant 
with |anti) <-> |11) causing a large error during the operation of the CNOT 
gate. While it is desirable to make J{t) large, it must be kept comfortably be- 
low 816.65 as near this level crossing the time required to adiabatically increase 
J(t) increases significantly. 
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Step 3 (the decreasing of A\) could be performed without degrading the 
overall fidelity of the gate in a time of less than a micro-second with a linear 
pulse profile. 

The above steps were simulated using an adaptive Runge-Kutta routine to 
solve the density matrix form of the Schrodinger equation in the computational 
basis \n\e\n-ie-i). 

p(t) = ^[H(t),p(t)]. (8) 
The times used for each stage are as follows 



stage 


duration (fis) 


2 


9.0000 


3 


0.1400 


4 


7.5989 


5 


9.0000 


6 


0.1400 



Note that the precision of the duration of stage 4 is required as the oscillating 
field -Bac induces the states |11) and |symm) to swap smoothly back and forth. 
The duration 7.5989/is is the time required for one swap. The other step times 
were obtained by first setting them to arbitrary values (~5/zs) and increasing 
them until the gate fidelity ceased to increase. The step times were then de- 
creased one by one until the fidelity started to decrease. As such, the above 
times are the minimum time in which the maximum fidelity can be achieved. 
This maximum fidelity was found to be 5 x 10 -5 for all computation basis states. 



4 Intrinsic dephasing and fidelity 

In this paper, dephasing is modelled as exponential decay of the off diagonal 
components of the density matrix. While a large variety of dephasing models 
exist |25l 1261 1771 1281 1231 130| . this approximation is consistent with the first order 
nature of the spin Hamiltonian. The donor electrons and phosphorous nuclei 
are assumed to dephase at independent rates. With the inclusion of dephasing 
terms (Eq. |HJ) becomes 

P = ^\H,R\ 

-r e [<,[<,p]]-r e [<,[<, P ]] 
-r n [<j<,p]]-r n [<j< )P ]]. (9) 

To understand the effect of each double commutator, it is instructive to consider 
the following simple mathematical example : 

M = -T[o*,[o*,M\] 
raw ™12 \ ( — AYmi2 



™2i m 2 2 J \ -4rr7i 2 i 
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( m u (t) m 12 (t) \ _ ( mn(0) m 12 (0)e~ 4r * \ 

^ m 21 (t) m 22 (t) ) ~ \ m 21 (0)e- 4r * m 22 (0) J ' liUj 

Thus each double commutator in (Eq.[!5J) exponentially decays it's associated off 
diagonal elements with a characteristic time r e = l/4r e or t„ = l/4r„. 

For each initial state |00), |01), 1 10} and |11) (Eq.|5J) was solved for a range 
of values of r e and r„ using the pulse profiles described in section [3] allowing 
a contour plot of the gate error versus r e and r„ to be constructed (Fig. |3J). 
Note that each contour is a double line as each run of the simulation required 
considerable computational time and the data available does not allow finer 
delineation of exactly where each contour is. The worst case error of all input 
states as a function of r e and r„ is shown in (Fig. |BJ. 



5 Conclusion 

At the time of writing, to the authors' knowledge the only experimental mea- 
surement of the 31 P in 28 Si dephasing times in is | 20| in which the donor electron 
dephasing time T 2 was measured at T=1.4K, B Z =0.3T to be ~0.5ms. The 28 Si 
sample contained (0.12 ± 0.08)% 29 Si (T^|. Note that a dopant concentration of 
~10 16 cm -3 was used implying a donor separation of ~50nm. If this T 2 is used 
for r e and if r„ is assumed to be a couple of orders of magnitude larger as in 
the case of the relaxation times, then from (Fig. EJ) the overall error probability 
would be just under 10 -3 

Theoretical calculation of r e and r n has been performed in for a 2D array 
of P donors spaced lOnm apart in pure 28 Si yielding r e = 2/xs and r„ = 10s. 
Such a short r e would imply an unacceptable error probability of about 10%. 
However, this same paper also contains similar calculations for natural silicon 
(4.7% 29 Si) with r e quoted as 200/xs which leads to an overall error probability 
just over 10~ 3 . The suppression of decoherence in this case arises from line 
broadening due to the presence of 29 Si nuclei. In the case of the Kane quantum 
computer, similar suppression can be achieved by biasing the A-electrodes such 
that nearby qubits have different spin-flip energies. Further investigation of this 
point is required. 

Though 10~ 3 is a large error probability, numerical simulations by Zalka sug- 
gest this may be tolerable |15| . Work is in progress on simulations to determine 
an acceptable error rate for the Kane architecture. 
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Figure 5: Probability of error during a CNOT operation as a function of r e and 
r n for input state a.) |00), b.) |01), c.) |10) and d.) 1 11) . The first qubit is the 
control. 
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Figure 6: The worst case probability of error during a CNOT operation as a 
function of r e and t„ for all input states. 
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